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We perform a diagrammatic calculation of the leading two- 
loop QCD corrections to the masses of the neutral CP-even 
Higgs bosons in the Minimal Supersymmetric Standard Model 
(MSSM). The results are valid for arbitrary values of the pa- 
rameters of the Higgs and scalar top sector of the MSSM. 
The two-loop corrections are found to reduce the mass of the 
lightest Higgs boson considerably compared to its one-loop 
value. The numerical results are analyzed in the GUT fa- 
vored regions of small and large tan (3 . Their impact on a 
precise prediction for the mass of the lightest Higgs boson is 
briefly discussed. 



Supersymmetric theories (SUSY) ffl are the best mo- 
tivated extensions of the Standard Model (SM) of the 
electroweak and strong interactions. They provide an 
elegant way to break the electroweak symmetry and to 
stabilize the huge hierarchy between the GUT and the 
Fermi scales, and allow for a consistent unification of 
the gauge coupling constants as well as a natural so- 
lution of the Dark Matter problem; for recent reviews 
see Ref. |J. The MSSM predicts the existence of scalar 
partners JltIr to each SM chiral fermion, and spin-1/2 
partners to the gauge bosons and to the scalar Higgs 
bosons. So far, the direct search for SUSY particles has 
not been successful. One can only set lower bounds of 
0(100) GeV on their masses §. 

A particularly stringent test of the MSSM is the search 
for the lightest Higgs boson. At tree level its mass, rrih, is 
predicted to be lower than the one of the Z boson. How- 
ever, the one-loop corrections are known to be huge H,|). 
As an impact, rrih > Mz is possible, and an upper bound 
of approximately 150 GeV is obtained. Hence a two-loop 
calculation is inevitable for a precise prediction of the 
mass of the lightest Higgs boson. This is particularly im- 
portant in view of the search for this particle at LEP2, 
where a precise knowledge of rrih in terms of the rele- 
vant SUSY parameters is crucial in order to determine 
the discovery (and of course also the exclusion) potential 
of LEP2. 

Up to now there exist renormalization group improve- 
ments of the one-loop result by including the two-loop 
leading logarithmic contributions [^[j8), and a diagram- 
matic calculation of the dominant two-loop contributions 
in the limiting case of vanishing i-mixing and infinitely 
large Ma and tan /3 M . These results indicate that the 



two-loop corrections considerably reduce the prediction 
for rrih- However, a diagrammatic two-loop calculation 
of the neutral mass spectrum going beyond the above- 
mentioned limiting case has been missing so far. Such 
a Feynman diagrammatic calculation is technically very 
involved, but it is of particular interest since it allows for 
general parameters of the MSSM Higgs sector and for vir- 
tual particle effects without restrictions on their masses 
and mixing. It is the purpose of this letter to investigate 
the leading QCD corrections to the masses of the neutral 
CP-even Higgs bosons and in particular to provide in this 
way a two-loop prediction of rrih for arbitrary values of 
the parameters of the Higgs and scalar top sector of the 
MSSM. 

Contrary to the SM, in the MSSM two Higgs doublets 
are needed. The Higgs potential is given by |l(J 
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where mi,m 2 ,mi2 are soft SUSY-breaking terms, g,g' 
are the SU(2) and U(l) gauge couplings, and £12 = —1. 

The doublet fields Hi and H 2 are decomposed in the 
following way: 
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The potential eq. (Q) can be described with the help of 
two independent parameters (besides g, g'): tan/? = 
v 2 /vi and M\ = — mf 2 (tan/3 + cot/3 ), where Ma is the 
mass of the CP-odd A boson. 

At tree level the mass matrix of the neutral CP-even 
Higgs bosons is given in the <pi — <p 2 basis in terms of Mz 
and Ma through 
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The tree-level mass predictions receive large correc- 
tions at one-loop order through terms proportional to 
GFrn\\r\{jai i mi^lm 2 ) These dominant one-loop con- 
tributions can be obtained by evaluating the contribution 
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of the t — i-sector to the </>i : 2 self-energies at zero exter- 
nal momentum from the Yukawa part of the theory (ne- 
glecting the gauge couplings). Accordingly, the one- loop 
corrected Higgs masses are derived by diagonalizing the 
mass matrix 



M 2 

Higgs 



S*j(0) r 
- £^ a (0) 



'0102 
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where the £ denote the Yukawa contributions of the t — t- 
sector to the renormalized one-loop 4>i,2 self-energies. By 
comparison with the full one-loop result || it has been 
shown that these contributions indeed contain the bulk 
of the one-loop corrections. They typically approximate 
the full one-loop result up to about 5 GeV. 

In order to derive the leading two-loop contributions to 
the masses of the neutral CP-even Higgs bosons we have 
evaluated the QCD corrections to eq. (|4|), which because 
of the large value of the strong coupling constant are 
expected to be the most sizable ones (see also Ref. ||). 
This requires the evaluation of the renormalized </)i,2 self- 
energies at the two-loop level. Typical Feynman dia- 
grams corresponding to the Yukawa contributions of the 
t — t-sector to the 0i j2 self-energies and tadpoles are 
shown in Fig. [I]. They have to be supplemented by the 
countcrtcrm insertions in the corresponding one-loop dia- 
grams. Fig. [j]a shows the pure scalar contributions to the 
Higgs self-energies. In Fig. [j]b the gluonic corrections are 
depicted, while Fig. shows the gluino-exchange contri- 
bution. In Fig. |l|d-f the tadpole contributions for these 
three types of corrections are given. 
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FIG. 1. Typical Feynman diagrams for the two- loop con- 
tribution to the Higgs-boson self-energies and tadpoles. 
H = 01,0a, A. 



The renormalization has been performed in the on- 
shell scheme. The counterterms in the Higgs sector are 
derived from the Higgs potential eqs. (Eh, (g) by expand- 
ing the counterterm contributions up to two- loop order. 
The renormalization conditions for the tadpole countert- 
erms have been chosen in such a way that they cancel the 
tadpole contributions in one- and two-loop order. The 
renormalization in the t — t-sector has been performed 
in the same way as in Ref. |fl]]. For the present cal- 
culation the one- loop counterterms Sm t , Sm^ , Sm^ 2 for 
the top-quark and scalar top-quark masses and SOi for 
the mixing angle contribute, which enter via the subloop 
renormalization. The appearance of the t mixing angle 
Ql reflects the fact that the current eigenstates, t^ and 
tR, mix to give the mass eigenstates t\ and £2- Since 
the non-diagonal entry in the scalar quark mass matrix 
is proportional to the quark mass the mixing is particu- 
larly important in the case of the third generation scalar 
quarks. The mixing-angle counterterm S9f is chosen such 
that there is no mixing between t\ and ti when t\ is on- 
shell. The numerical result, however, is insensitive to 
this choice of the renormalization point. The one-loop 
counterterms for /i and tan (3 , Sfj, and 5 tan (3 , do not 
contribute since they are independent of a s . 

The renormalized self-energies have the following 
structure: 



E,(0) = E«(0) + £( 2 >(0) - SVj^ - SV} 2 \ 



(5) 



where s = 4>i,4>2,4>i4>2- £» and Eg 2 -* denote the un- 
renormalized self-energies at the one- and two-loop level, 
and SVj 1 ^ and SVj 2 ^ are the one- and two-loop countert- 
erms derived from the Higgs potential. The counterterms 
read: 
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Ta', where T® denotes the tadpole con- 
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tribution, St® is the corresponding counterterm, and 
5M 2 ® = E®(0) (z = 1,2). 

In deriving our results we have made strong use of 
computer-algebra tools. The package FeynArts jl^] (in 
which the relevant part of the MSSM has been imple- 
mented) has been applied to generate the Feynman am- 
plitudes and the counterterm contributions. For evalu- 
ating the amplitudes the package TwoCalc |l5| ] has been 
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used. The calculations have been performed using Di- 
mensional Reduction (DRED) (m), which is necessary in 
order to preserve the relevant SUSY relations. Naive ap- 
plication (without an appropriate shift in the couplings) 
of Dimensional Regularization (DREG) |lq| , on the other 
hand, does not lead to a finite result. The same observa- 
tion has also been made in Ref. |)). 

The contributions of the scalar, the gluon-, and the 
gluino-exchange diagrams in Fig. [I] together with the cor- 
responding counterterm contributions are not separately 
finite (as it was the case in the calculation of Ref. [pd|), 
but have to be combined in order to obtain a finite re- 
sult. Our results for the two- loop self-energies are 
given in terms of the SUSY parameters tan/3 , Ma, /i, 
m% , TO£ 2 , Of, and m g . In the general case the results are 
by far too lengthy to be given here explicitly. In the spe- 
cial case of vanishing mixing in the i-sector, /1 = 0, and 
a relatively compact expression can be 
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derived. It is given by 
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rj being the solutions of to|t + m 2 /r = m 2 + to 2 — p 2 
(j = 1,2). Eq. (|l(]| ) approximates the complete numerical 
result for vanishing mixing (for arbitrary [i and mj ^ 
TOj 2 ) up to about 2% accuracy. 

Inserting the one- loop and two- loop ^2 self-energies 
into eq. (0), the predictions for the masses of the neu- 
tral CP-even Higgs bosons are derived by diagonaliz- 
ing the two-loop mass matrix. For the numerical eval- 
uation we have chosen two values for tan/3 which are 
favored by SUSY-GUT scenarios tan/3 = 1.6 for 

the SU(5) scenario and tan/3 = 40 for the 50(10) sce- 
nario. Other parameters are Mz — 91.187 GeV, M\y = 
80.375 GeV, G F = 1.16639 10~ 5 GeV" 2 , a s = 0.1095, 
and TOt = 175 GeV. For the figures below we have fur- 
thermore chosen fi = -200 GeV, M A = 500 GeV, and 
rn g = 500 GeV as typical values. The scalar top masses 
and the mixing angle are derived from the parameters 

LR 



Mi and M t of the t mass matrix (our conven- 



tions are the same as in Ref. |ll[) 
we have chosen m r ~, = Mj = Mf . 



In the figures below 



The plot in Fig. | shows to/j as a function of M^ R j r 



where m g is fixed to 500 GeV 



A minimum is reached 
for M^ R — GeV which we refer to as 'no mixing'. A 
maximum in the two-loop result for to/j is reached for 
about M t Lfl /TOg « 2 in the tan/3 =1.6 scenario as well 
as in the tan (3—40 scenario. This case we refer to 
as 'maximal mixing'. Note that the maximum is shifted 
compared to its one-loop value of about M/^/to, w 2.4. 
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FIG. 2. One- and two-loop results for mh as a function of 
Mt R juiq for two values of tan/3 . 
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In Fig. H the low-tan scenario with tan/3 = 1.6 is 
analyzed. The tree-level, the one-loop and the two-loop 
results for rrih are shown as a function of rriq for no mixing 
and maximal mixing. For both cases the one-loop result 
is in general considerably reduced. For the no-mixing 
case the difference between the one-loop and two-loop 
result amounts up to about 18 GeV for rriq = 1 TcV. In 
the maximal-mixing case the reduction of the one-loop 
result is about 10 GeV for m q = 260 GeV (for smaller rriq 
one gets unphysical or experimentally excluded t-masses) 
and about 25 GeV for m 3 = 1 TeV. 
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FIG. 3. The mass of the lightest Higgs boson for 
tan/3 = 1.6. The tree-, the one- and the two-loop results 
for mh are shown as a function of rriq for the no-mixing and 
the maximal-mixing case. 

The variation of this result with m g is of the order of 
few GeV. Varying tan/3 around the value tan/3 =1.6 
leads to a relatively large effect in rrih- Higher values 
for rrih are obtained for larger tan/? . A more detailed 
analysis of the dependence of our results on the differ- 
ent SUSY parameters will be presented in a forthcoming 
publication. 

In Fig. U the high-tan f3 scenario with tan (3 — 40 
is analyzed. Again the tree-level, the one-loop and the 
two-loop results for rrih are shown as a function of rriq for 
minimal and maximal mixing. As in the case of low tan /?, 
the one- loop result is in general considerably reduced. 
For no mixing the difference between the one-loop and 
two-loop result reaches about 14 GeV for rriq = 1 TeV. 
In the maximal-mixing case the reduction of the one-loop 



result amounts to about 7 GeV for nig 



260 GeV and 



about 22 GeV for m; 



1 TeV. The reduction of the 



one-loop result is slightly smaller than for tan/3 = 1.6. 
This can be understood from the result for E^, 2 (0) given 
as a special case in eq. (M. In this case (3 appears only 
in the prefactor as 1/ sin (3 and one thus gets a bigger 



reduction of rrih for smaller tan/3. The variation of the 
result shown in Fig. ^ with nig is again of the order of 
few GeV. The effect of varying tan (3 around tan (3 = 40 
is marginal. 
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FIG. 4. The mass of the lightest Higgs boson for 
tan/3 = 40. The tree-, the one- and the two-loop results 
for rrih are shown as a function of rriq for the no-mixing and 
the maximal-mixing case. 

We have compared our results with the results ob- 
tained in Ref. |j] in the case of no i-mixing and Ma — > 
oo, tan /3 — > oo and have checked analytically that in the 
limiting case = = m g 3> m t in eq. (Mh we re- 
cover the corresponding formula given in Ref. 

Supplementing our results for the leading 0(aa s ) cor- 
rections with the leading higher-order Yukawa term of 
C(a 2 mf) given in Ref. M leads to an increase in the 
prediction of rrih of up to about 3 GeV. A similar shift 
towards higher values of rrih emerges if at the two-loop 
level the running top-quark mass, mt(rrit) = 166.5 GeV, 
is used instead of the pole mass, m t = 175 GeV, thus 
taking into account leading higher-order effects beyond 
the two-loop level. We have compared our results with 
the results obtained by two-loop renormalization group 
calculations given in Refs. §,§]] We find good agree- 
ment for the case of no i-mixing, while for larger i-mixing 
sizable deviations exceeding 5 GeV occur. In partic- 
ular, the value of M^ R /rriq yielding the maximal rrih 
is shifted from M^ R jrriq w 2.4 in the one- loop case to 
Mi R / rriq ss 2 when our diagrammatic two-loop results 
are included (see Fig. ||). In the results based on renor- 



1 The results of Ref. 
2 GeV with each other. 



and Ref. M agree within about 
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malization group methods ||]8), on the other hand, the 
maximal value of is obtained for M^ R /rriq w 2.4, i.e. 
at the same value as for the one- loop result. 

In summary, we have diagrammatically calculated the 
leading 0(aa s ) corrections to the masses of the neutral 
CP-even Higgs bosons in the MSSM. We have applied 
the on-shell scheme and have imposed no restrictions on 
the parameters of the Higgs and scalar top sector of the 
model. The two-loop correction leads to a considerable 
reduction of the prediction for the mass of the lightest 
Higgs boson compared to the one-loop value. The reduc- 
tion turns out to be particularly important for low values 
of tan [3 . The leading two-loop contributions presented 
here can directly be combined with the complete one- 
loop results in the on-shell scheme ||. A discussion of 
the corresponding results will be given in a forthcoming 
paper, where also a more detailed comparison with the 
results based on renormalization group methods will be 
pursued. 
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